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ABSTRACT 

From a deep multi-epoch Chandra observation of the elliptical galaxy NGC 3379 we report the spectral 
properties of eight luminous LMXBs (Lx> 1-2 x 10 38 erg s" 1 ). We also present a set of spectral simulations, 
produced to aid the interpretation of low-count single-component spectral modeling. These simulations demon- 
strate that it is possible to infer the spectral states of X-ray binaries from these simple models and thereby 
constrain the properties of the source. 

Of the eight LMXBs studied, three reside within globular clusters, and one is a confirmed field source. Due 
to the nature of the luminosity cut all sources are either neutron star binaries emitting at or above the Eddington 
luminosity or black hole binaries. The spectra from these sources are well described by single-component 
models, with parameters consistent with Galactic LMXB observations, where hard-state sources have a range 



2 



Brassington et al. 



in photon index of 1 .5-1 .9 and thermally dominant sources have inner disc temperatures between ~ 0.7 - 1 .55 
keV. 

The large variability observed in the brightest globular cluster source (Xx> 4 x 10 38 erg s" 1 ) suggests the 
presence of a black hole binary. At its most luminous this source is observed in a thermally dominant state 
with fcTi„=1.5 keV, consistent with a black hole mass of ~4M . This observation provides further evidence 
that globular clusters are able to retain such massive binaries. We also observed a source transitioning from 
a bright state (Lx^ 1 x 10 39 erg s -1 ), with prominent thermal and non-thermal components, to a less luminous 
hard state (Lx=3.8x 10 38 erg s , T =1.85). In its high flux emission this source exhibits a cool-disc component 
of ^0.14 keV, similar to spectra observed in some ultraluminous X-ray sources. Such a similarity indicates a 
possible link between 'normal' stellar mass black holes in a high accretion state and ULXs. 
Subject headings: galaxies: individual (NGC 3379) — X-rays: galaxies — X-ray: binaries 



1. INTRODUCTION 

The discovery with Chandra of several low-mass X-ray 
binary (LMXB) populations in early-type galaxies, and the 
associations of these LMXBs with either Globular Clusters 
(GCs) or the stellar field, have provided new impetus to the 
study of the formation and evolution of LMXBs in GCs and 
to the possible relation of field LMXBs to the GC population 
(Grindlay & Hertz 1985; Verbunt & van den Heuvel 1995; 
see Fabbiano 2006 and refs. therein). Given the characteris- 
tics of these data, most of this work has been based on pop- 
ulation studies (e.g., Kim et al. 2006: Kundu et al. 2007; 
Sivakoff et al 2007; Voss & Gilfanov 2007). However, in the 
few cases of detection of luminous sources in deep enough 
observations, detailed spectral and variability studies can be 
pursued, to provide more direct constraints on the nature of 
the X-ray sources. One such example is the recent discovery 
of a variable luminous GC source in NGC 4472, with tempo- 
ral and spectral characteristics supporting a stellar BH binary 
(Maccarone et al 2007; Shih et al 2008). 

In this paper we report the results of the spectral analysis of 
eight luminous sources (Lx> 1 -2 x 10 38 erg s -1 in the 0.3-8.0 
keV band), detected in the nearby elliptical galaxy NGC 3379 
(in the poor group Leo, D=10.6 Mpc, Tonry et al. 2001 ') 
with Chandra ACIS-S (Weisskopf et al 2000). These sources 
were observed at five different epochs, as part of a monitoring 
campaign with Chandra (PI: Fabbiano) providing the rare op- 
portunity of long-term spectral monitoring of LMXBs in an 
elliptical galaxy. These sources are part of the sample of 132 
sources detected in NGC 3379 from these observations with 
luminosities greater than a few 10 36 erg s" 1 (Brassington et 
al. 2008 - hereafter B08). 

In section|2]we describe the observations and the properties 
of the sources under study from the B08 catalog. In section 
[3] we describe our spectral analysis and report the results. In 
section [4] we present spectral simulations that we have per- 
formed to aid in the interpretation of our results and in section 
|5] we discuss the spectral analysis and compare the results to 
our simulations. Our conclusions are summarized in section 
i 

2. SOURCE SELECTION 

From the catalog of B08 we have selected a sub-set of the 
brightest sources (Lx> 1-2 x 10 38 erg s" 1 ) within the D25 el- 
lipse of the galaxy to perform detailed spectral and temporal 
analysis. This luminosity cut results in the selection of 12 
sources, we further exclude the nuclear source (S81) and the 

Electronic address: nbrassington@head.cfa.harvard.edu 

1 Note that Jensen et al. (2003) report a distance of 9.82 Mpc. How- 
ever, the distance of Tonry et al. (2001) is adopted here for consistency with 
Brassington et al. (2008). This choice does not affect our conclusions in any 
way. 



confirmed ULX (S70; Fabbiano et al. 2006), both of which 
will be the subject of forthcoming papers. We also exclude 
confused sources in the central region (S75 & S82), where 
reliable spectral extraction cannot be carried out. This re- 
sults in a final selection of eight sources. Of these, six lie 
within the HST WFPC2 field of view of NGC 3379, allow- 
ing firm constraints to be placed on their optical counterparts 
(see B08). Three of these sources (S41, S42 & S67 of B08) 
are identified with GCs and one has been confirmed as a field 
LMXB (S86). The other two sources are found within the 
region of confusion from the HST observation (S74 & S77) 
and are therefore unclassified. The two remaining sources 
are external to the optical FOV (S102 & S103) and are there- 
fore also unclassified. At a flux of 8.9xl0~ 15 erg s" 1 cm -2 
(which gives L x =1.2x 10 38 erg s" 1 at D=10.6 Mpc), from the 
ChaMP + CDF logN-logS relation of Kim et al. (2007), 
less than 1 source is expected to be a background AGN over 
the area encompassing the eight sources. The Lx> 1.2 x 10 38 
erg s -1 luminosity cut ensures that the selected sources have 
a minimum of ^400 counts from the co-added observation 
and 100 counts in at least one individual observation, allow- 
ing meaningful spectral analysis to be carried out. By nature 
of the selection criteria these sources are either NS binaries 
emitting near or above the Eddington luminosity, or are black 
hole binaries. 

From the observations of NGC 3379 the parameters de- 
rived in B08 indicate that BH-LMXBs are present in the X-ray 
source population of this galaxy, both in the field and in GCs. 
BH-LMXBs in the field are a well-established class of X-ray 
binaries in the Galaxy (see the review of Remillard & McClin- 
tock 2006; from hereon in RM06). Although BH-LMXBs are 
not expected to be common or luminous in GCs (Kalogera 
et al. 2004; see Ivanova et al. 2010 for a recent discussion 
of the formation of BHs within GCs), there is recent con- 
vincing evidence of at least two such luminous sources. The 
first source, a ULX detected within the elliptical Virgo galaxy 
NGC 4472, exhibits strong X-ray variability (Maccarone et 
al. 2007, Shih et al. 2008) providing unambiguous evidence 
for a BH-GC. Subsequent optical spectroscopy of the GC has 
revealed very broad [O III] emission lines, leading to the sug- 
gestion that this source is a stellar mass black hole accreting 
above or near the Eddington limit (Zepf et al. 2008). The 
second source, a ULX in NGC 1399 (Irwin et al. 2010), also 
exhibits strong [O III] emission lines (although these are less 
broad than the NGC 4472 source) and has little or no hydro- 
gen emission. This source could also be a stellar mass black 
hole accreting at or near the Eddington limit, although Irwin 
et al. (2010) suggest that tidal disruption of a white dwarf 
from an intermediate-mass black hole (IMBH: 100-10 4 M ) 
could explain the optical emission of this GC. From this in- 
terpretation a minimum black hole mass of 1000 M Q is im- 
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plied. If sources S41, S42 & S67 are luminous BH binaries, 
the estimated a-posterior probability of a BH-GC association 
is ~4%, given that 70 GCs are found in NGC 3379 within 
the joint Chandra/Hubble field of view (see B08 and refs. 
therein). Given that only 10 GCs are associated with X-ray 
sources in NGC 3379 (for L x > a few 10 36 erg s" 1 ; B08 and 
Kim et al. 2009), these luminous BH-LMXBs may constitute 
a significant fraction of the GC-LMXB population. 

The eight sources presented here are shown in Figure Q] 
where an adaptively smoothed 0.3-8.0 keV X-ray image is 
presented, with the source regions and the HSJ FOV overlaid. 
In Table Q] the log of the Chandra observations is reported, 
along with the background-subtracted net source counts in 
the 0.3-8.0 keV band, for each source in each observation 
(as reported in B08). In this table the optical correlations 
are also identified, where GC indicates a confirmed globular 
cluster counterpart, F a field sources and - indicates sources 
with insufficient optical data to confirm a counterpart (or lack 
thereof). 

Figure |2] (adapted from B08) shows the light curves of the 
eight sources from the Chandra observations with panels in- 
dicating X-ray luminosity, the hardness ratio (HR) and color 
values of the source in each observation. In the top panel 
Lx indicates the X-ray luminosity in the 0.3-8.0 keV band, 
for an assumed power law spectrum with T=1.7 and Galactic 
line-of-sight 2 , where the horizontal dashed line indicates the 
Lx derived from the co-added observation. Because the same 
canonical model has been applied to each source this value 
represents a scaling of the extracted source count-rate, where 
both the spatial and the temporal quantum efficiency varia- 
tions have been accounted for by generating an energy conver- 
sion factor for each sources in each observation (see §2.1 in 
B08 for more details). In the second panel, the HR variation 
is shown where HR=(Hc-Sc)/(Hc+Sc), with Sc and He net 
counts in the 0.5-2.0 keV and 2.0-8.0 keV bands respectively. 
In the third and fourth panels the X-ray colors C21 and C32 
are presented, these colors are defined as C21 = log(Si/S2) 
and C32 = log(S2/H), where Si, S2 and H are the net counts 
respectively in the energy bands of 0.3-0.9 keV, 0.9-2.5 keV 
and 2.5-8.0 keV. As can be seen from these figures, the eight 
sources are persistent, in that they were detected with compa- 
rable luminosity over a five year span, although the majority 
exhibits some variability in luminosity and/or spectral prop- 
erties. 

In Table|2]the long- and short-term variability of the sources 
are summarized, where the short-term variability is defined 
for each pointing, by means of the Kolmogorov-Smirnov 
test (K-S test). Here V indicates variable sources (with val- 
ues >99% confidence), P indicates possible variable sources 
(with variability values >90% confidence) and the - indicates 
non-variable sources. The long-term variability from all five 
observations was defined by the chi-squared test where V in- 
dicates variability, and N indicates a non-variable source. In 
the last column, Sign presents the significance in change in 
luminosity between the highest and lowest flux observations. 
All of these classification are fully described in B08, §2.4. All 
the sources, with the exception of S67, have been determined 
to exhibit long-term variability. 



3. SPECTRAL ANALYSIS AND RESULTS 

2 Ny l =2.1%X\O l0 cm-' 2 calculated with the tool COLDEN: 
http://cxc.harvard.edu/toolkit/colden.jsp 



The data were processed and calibrated as described in B08. 
Analysis was performed with the CXC CIAO software suite 
(v3.4) 3 with CALDB version 3.5 and HEASOFT (v5.3.1). 
Spectra were extracted for every pointing for each of the eight 
sources, using the CIAO tool psextract, where source regions 
were defined to be circular regions with radii as in B08; back- 
ground counts were extracted from surrounding annuli, with 
outer radii 2-3 times larger, depending on the presence of 
nearby sources. In cases where the source region was found 
to overlap with a nearby sources, the extractions region ra- 
dius was reduced and the area of the overlapping source was 
excluded. The minimum defined extraction radius was 1.5" 
on-axis, and 3" off-axis (maximum off-axis angle ~1.5'), en- 
suring that the source regions enclosed >93% of the encircled 
energy in all cases. 

The source spectra were fitted in XSPEC (vll.3), where 
the data were restricted to 0.3-8.0 keV, as energies below this 
have calibration uncertainties, and the spectra presented here 
do not have significant source flux above 8.0 keV, where the 
data are highly contaminated by cosmic rays. The spectra 
were fitted to two models that describe the properties of X- 
ray binary spectra well (see e.g. review by RM06): the multi- 
color disc blackbody (DISKBB in XSPEC; hereon in referred 
to as DBB), and the power law (PO) models, along with the 
wabs photoelectric absorption model 4 . All parameters were 
allowed to vary freely, with the exception of the absorption 
column A^h, which in cases where the best-fit value was be- 
low that of the Galactic absorption, was frozen to the Galactic 
value of 2.78 x 10 20 cm" 2 . In instances where there were suffi- 
cient counts to bin the data to at least 20 counts per bin (allow- 
ing Gaussian error approximation to be used), the minimum 
X 2 method was used to fit the data. Where there were too 
few counts, the Cash statistic (Cash 1979) was used in prefer- 
ence to x 2 , where the spectra were extracted without binning. 
However, this statistical method has the disadvantage that it 
does not provide a goodness-of-fit measure like x 2 . 

Spectra for each source were first fitted separately for each 
pointing. These individual best-fit models, alongside the Lx 
and HR behavior shown in Figure [2] were then used to guide 
the joint fits for each source. In Table [3] these co-added fits 
are summarized, where column (1) gives the source number, 
column (2) the observations used in each fit, column (3) the 
net source counts, column (4) indicates which spectral model 
was used, column (5) the fit statistic (x 2 and number of de- 
grees of freedom v, or the C statistics - indicated by C), col- 
umn (6) the null hypothesis probability (or goodness when 
using the Cash statistic). Columns (7), (8) & (9) present the 
best-fit values of the fit parameters with 1 a errors for 1 inter- 
esting parameter: A^h (2.8Fr denotes that the value was frozen 
at Galactic column density), T for the PO model and fcT; n , the 
temperature of the innermost stable orbit of an accretion disc, 
for the DBB. Column (10) indicates the intrinsic value of Lx 
(calculated as a weighted average of the individual luminosi- 
ties from the joint fit) and column (11) the source luminosity 
range. This range was taken from the (ler) lower-bound value 
of the lowest flux observation included in the joint spectral 
fit, to the (ler) upper-bound value of the highest flux obser- 
vation included in the joint fit. When only one spectrum was 
modeled, the range indicates the lower- and upper-bound Lx 
values from that single observation. 

3 http://asc.harvard.edu/ciao 

4 Given the statistics of our data and the ACIS resolution, use of other 
absorption models does not affect our results and conclusions (see j)4.3t . 
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3.1. Globular Cluster Sources 

Of the sub-sample of sources presented in this paper, three 
have been determined to be GC-LMXBs. Through inspec- 
tion of their L x , HR and color values presented in Figure |2l 
and the individual spectral fits of each source, two LMXBs 
(S41 & S67) have been determined to exhibit little flux and 
no spectral variability, whilst the remaining source (S42) has 
been found to exhibit spectral changes between each observa- 
tion, with a variance in flux significance of 6a. In the case of 
the non-varying LMXBs, the five spectra for each source were 
jointly fitted in XSPEC, where all parameters, apart from the 
normalizations, were linked. For both sources the DBB and 
PO models were fitted to the spectra, these best-fits are sum- 
marized in Table [3] In cases where the DBB model was sta- 
tistically rejected the PO model only is presented. 

For the variable source, S42, the most luminous GC-LMXB 
in our sample, the individual fits and Lx and HR values from 
B08 indicate that the source has been observed in three differ- 
ent spectral states. The spectra have therefore been grouped 
to reflect these differences. The first grouping contains spec- 
tra with lower Lx and soft HR value (obs 1 &3), the second, 
spectra with higher Lx and hard HR value (obs 2&4) and the 
remaining spectrum indicates an intermediate state (obs 5). 
Both the DBB and PO models were applied to these three 
joint fits and are summarized in Table[3] The 90% confidence 
contours for two interesting parameters are presented in Fig- 
ure [3] where the left-hand panel shows the PO model in the 
three spectral groups and right-hand panel the DBB fits. In 
both panels Galactic Nh is indicated by the vertical solid line. 

3.2. The Field Source 

Only one source presented in this paper has been confirmed 
as a field-LMXB, S86. It has been determined that this source 
not only exhibits long-term variable behavior in both its flux 
and spectra, but that it also displays short-term flux variability 
(defined as a possible variable source with the K-S test: Table 
I2J) during the second pointing. This short-term variability has 
been further investigated through examining the light-curve 
of the source during observation 2, which was created with 
the CIAO tool dmextract. To provide well constrained val- 
ues of count-rate, allowing variability to be identified, binning 
of 10000 s was used when extracting counts. This resulting 
light-curve is shown in Figure |4j where it can be seen that the 
count-rate drops from ^0.008 cnt s" 1 to ^0.006 cnt s" 1 part- 
way through the observation. Spectra were extracted from 
these two different count-rate epochs, to identify if a change 
in spectral shape accompanied this change in flux, however, 
it was determined that, within errors, the best-fit values from 
the two spectra were identical. 

Once it had been determined that the short-term variabil- 
ity of S86 did not correspond to any spectral variation, the 
long-term spectral variability was investigated. Once again, 
guided by the individual spectral fits and the X-ray luminos- 
ity and HR values, the spectra were separated into two groups; 
a low Lx hard HR state (1,3&5) and a high Lx softer HR state 
(2&4). These best-fit values are summarized in Tableland 
their confidence contours are presented in Figure[5] 

3.3. Unclassified Sources 

The remaining four sources in this paper have insufficient 
optical data to confirm them as field or GC sources, either 
residing with the central 5" of the galaxy, where source con- 
fusion means GC sources cannot be easily identified, or are 



detected outside of the HST FOV. Of these sources, two (S77 
& S103) have been determined to have no spectral variability, 
despite both of them being being classified as having flux vari- 
ability. Consequently all five observations were fitted jointly 
for each source, following the methods described in 33. II for 
S41 & S67. 

Of the two remaining LMXBs (S74 & S102), both have 
been determined to have variability and require separate spec- 
tral modeling. S74 exhibits stable long-term flux and spec- 
tra over the first four observations but in observation five a 
four-fold increase in flux (with a significance in flux change 
of 16a) and a softening in HR is observed. The joint fits 
have therefore been separated into obs 1,2,3&4, and obser- 
vation 5, to reflect these variations. When modeling the spec- 
trum from observation 5 it was found that there was a drop 
in count-rate in the binned data at ~ 1 .15 keV and from ad- 
ditional analysis it was determined that this feature was still 
present in the unbinned data. To further investigate this, addi- 
tional components were included in the spectral model in an 
attempt to describe this feature. However, no physical mod- 
els described the data well. Background subtraction, edge ef- 
fects and calibration uncertainties were also considered as an 
explanation of the low energy bin but after further investiga- 
tion all were rejected. After ruling these possibilities out, the 
most likely remaining explanation of this 'feature' is that it is 
a consequence of statistics, and not a physical feature of the 
source. The PO and DBB models were fitted to the whole of 
the 0.3-8.0 keV spectrum and these best-fit model parameters 
are presented in Table [3] In subsequent fits the energy range 
of 1.1—1.2 keV was excluded to determine the effect of the 
low-count bin. From this it was shown that the best-fit values 
are consistent with those from the full spectrum fit, although 
the goodness-of-fit values are greatly improved by excluding 
this energy range (PO y^lv from 1 .69 to 1 .30; DBB j^lv from 
1.62 to 1.26). Some short-term (possible) variability was also 
observed in observation 1 of S74, although there are too few 
counts in this pointing to be able to further investigate this. 

Through inspection of not only the long-term light-curve 
and HR values of SI 02 (Figure [2]) but also its individual fits, 
there is a suggestion of a gradual and continual change in flux. 
Initially this is an increase in flux from the archival observa- 
tion to the more recent pointings, followed by a steady de- 
cline, which is coupled with a suggestion of hardening in HR 
values. Therefore the joint fits were divided into three groups. 
Observation 1, the softest HR value, observations 2&3, high 
Lx, and observations 4&5, lower Lx- As with the previous 
sources, single-component PO and DBB models were applied 
to the spectra. 

However, for group 2, when the source count-rate is at its 
highest, the best-fit values from the PO model did not provide 
physically realistic parameters where the best-fit value of Nu 
was significantly below Galactic Nh, therefore a composite 
DBB+PO model was applied to the data. 

The best-fit spectra of these four sources are all summa- 
rized in Table [3] and the confidence contours of S74 are pre- 
sented in Figure [6] The confidence contours from the single- 
component PO model of the three spectral groupings from 
SI 02 are presented Figure [7] and the contours from the two- 
component model of observations 2 & 3 are shown in Figure 
i 

4. SIMULATIONS OF LMXB SPECTRA 

From studies of Galactic binaries it has been observed that 
the energy spectra of LMXBs often exhibit a composite shape 
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of both thermal and non-thermal emission components. In 
these sources the thermal component is well modeled by the 
multicolor disc blackbody model, which originates in the in- 
ner accretion disc, and the non-thermal component is well de- 
scribed by a power law component. This composite emission 
provides a wide range in X-ray properties that can be clas- 
sified into five distinct spectral states. In RM06 these states 
are defined to be; the quiescent state, a non-thermal hard state 
with very low emission (L x 10 3() 5 -10 33 5 erg s -1 ), the hard 
state (denoted hereon in as NT; non-thermal state, to avoid 
confusion with the high/soft state classification that can also 
be used to describe LMXB spectral states), where the non- 
thermal component dominates the spectra (>80% 5 ) and the 
spectral index is in the range 1.5< T < 2.1, and the thermally 
dominant state (TD), where the thermal emission arising from 
the disc contributes to more than 75% of the flux. In the TD 
state the disc temperature is typically in the range of 0.7-1.5 
keV and the faint power law component is steep (T >2.1). 
Sources have also been observed in the steep power law (SPL) 
state, where the sources are much brighter (Lx>0.2LEdd) with 
a significant flux contribution arising from the non-thermal 
component with steeper slope of >2.4, compared to ~1.7 in 
the hard state. The final spectral class encompasses emission 
not defined by the other states, where spectra appear to show 
sources between states, this is termed the intermediate state. 
It should be noted that this definition does not represent a sin- 
gle state but indicates that the source emission cannot be de- 
scribed by any of the four previous source classes. 

From Chandra and XMM-Newton observations of extra- 
galactic LMXBs there are typically too few counts to apply 
composite models to each spectrum, therefore single compo- 
nent models that provide an approximation of the spectra must 
be applied to the data. To understand how well these single- 
component models describe the X-ray emission of LMXBs, 
simulations have been performed with Sherpa (v3.4), using 
Chandra response files. In these simulations spectra have 
been produced from two-component models comprising both 
a MCD (XSDISKBB in Sherpa; the XSPEC DBB model used 
in section O and a PO component, along with the phabs pho- 
toelectric absorption model which has a value of Galactic ab- 
sorption. This model covers the 0.3-8.0 keV Chandra en- 
ergy range. These generated spectra have then been fitted to 
single-component models to characterize the typical best-fit 
values that are recovered. To broadly cover physical proper- 
ties that have been observed in the spectra of Galactic LMXBs 
(RM06) simulations were run over a large parameter space 
where the inner disc temperature of the DBB model ranged 
from 0.50 keV to 2.25 keV in steps of 0.25 keV and input Y 
slopes increased in steps of 0.2, from 1.5 to 2.5. Further, to in- 
vestigate how the composite spectral shape of sources in hard, 
thermally dominant and steep power law states were repre- 
sented by single-component fits, different flux ratios of 75% 
,60%, 50%, 40%, 25%, 10% and 0% were also included in the 
simulations, where the percentage indicates the flux contribu- 
tion from the thermal component. This was also expanded 
to include a sub-set of parameters at 90% and 100%, where 
spectra with T fixed at 1.7 and 2.5 were produced for the full 
disc temperature range. It was necessary to produce only a 
sub-set of these flux ratio values due to the computer inten- 

5 Where the full-band flux is 2-20 keV from RXTE spectra, not the 0.3-8.0 
keV band that has been used in the rest of this paper. Although, definitions 
from RM06 do provide a fairly accurate interpretation of spectral states ob- 
served in the Chandra bandpass. 



sive nature of these simulations. A value of 2.5 was selected 
for r based on typical LMXB photon index values of sources 
in a thermally dominant state. A photon index of 1.7 was 
selected based on the results from the lower flux ratio simu- 
lations, where this input parameter was found to provide me- 
dian values of A^h, T and k1 m for the single component best-fit 
parameters over the full photon index range. 

To replicate the quality of data that has been presented in 
this paper, the generated models were scaled to produce spec- 
tra of 1000 counts, with subsequent runs also covering a sub- 
set of 500 and 250 count data. After a spectrum was generated 
in Sherpa, using the fakeit 6 command, single-component PO 
and DBB models were fitted and the best-fit values recorded. 
Following this, the original two-component model was fitted 
to the data to ascertain how reliably the low-count spectrum 
could recover the input parameters. For each set of parame- 
ters 100 spectra were generated, providing a measure of the 
standard deviation of the best-fit values. 

4.1. Single Component Spectra 

4.1.1. Single PO Fits 

In Figure [9] a summary of the results of fitting the single- 
component PO model is presented. In the left-hand panel 
the best-fit value of A^h for each of the parameter models is 
shown and in the right-hand panel the difference (A) between 
input and best-fit photon index is presented. In both panels 
the best-fit values are plotted against increasing thermal flux 
ratios, with the solid diagonal lines indicating the best-fit val- 
ues for T=1.7, for each temperature step of the disc compo- 
nent (a sub-set of temperatures are shown for clarity in the 
figure). The shaded regions indicate the maximum and min- 
imum value of A^h over the full photon index range for each 
temperature step (where only photon index values of 1 .7 & 2.5 
have been simulated for flux ratio values of 90% and 100%). 
In the left-hand panel the horizontal dashed line indicates the 
value of Galactic A^h for NGC 3379. The standard deviation 
for the 1000 count spectra for the 75% and 25% flux ratios 
are indicated by the error bar in the top left corner of the 
plot, where the smaller error bar indicates the 25% flux ra- 
tio standard deviation (<r=1.82x 10 20 cm~ 2 ) and the larger er- 
ror bar the standard deviation from the 75% flux ratio model 
(o-=3.30xl0 20 cm- 2 ). 

From these figures it is clear that as the spectra become 
more disc-dominated the best-fit values from the single com- 
ponent PO model increasingly diverge from the input param- 
eters of the two-component model. Importantly, for all values 
of kT m and T, best-fit values of A^h are above Galactic (except- 
ing the spectra generated from the extreme models with input 
parameters T=2. 3-2.5 and kT m > 1 .50 keV) and even up to ~8 
times that of Galactic when the spectra are in a thermally- 
dominant state (i.e. the disc fraction contributes more than 
75% of the total flux). In fact, assuming a typical disc temper- 
ature of 1 keV, even with a disc contribution of only 40% the 
recovered A^h from the single-component PO fit is >3 times 
Galactic A^h (a value significantly above the Galactic value 
even after considering the standard deviation). Further, unsur- 
prisingly, the difference between the input and best-fit photon 
index (I ATI) increases as the ratio of disc flux to total flux 
increases, indicating that, as the non-thermal component be- 
comes less prominent within the composite fit, the reliability 
of the r recovered from the single component model is re- 
duced. 

6 http://asc. harvard. edu/sherpa3.4/threads/fakeit/ 
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A sub-set of these simulations were also generated for spec- 
tra with 500 and 250 counts, where flux ratios models of 
100%, 90%, 75%, 60%, 50%, 25% and 0%, were produced 
with an inner disc temperature of 1 keV and photon index 
of 1.7. In addition, for the higher flux ratio values >60%, 
models with T=2.5 were also used as input parameters for the 
simulations to represent thermally dominant states. The best- 
fit values of A^h with lower count spectra are compared with 
the full 1000 count simulations in Figure [TUJ where the 1000 
count spectral results are presented in an identical manner to 
the left-hand panel Figure [9] The best-fit A^h values of the 
500 and 250 count spectra are indicated by the blue and green 
lines respectively, with the shading indicating the best-fit val- 
ues from the T=2.5 simulations. As in Figure [9] the standard 
deviation of the 75% and 25% flux ratios are indicated by the 
color-coded error bars in the top left corner of the plot. From 
this figure it is clear that, despite the larger standard devia- 
tions in the lower-count models, the best-fit values of A^h are 
almost identical to those obtained with the 1000 count data. 
Even in the case of 250 counts, when a LMXB is in a 'typi- 
cal' thermal-dominant state (i.e. kT\ n =l keV, T=2.5 and disc 
flux >75%) the best-fit value of A^h is still more than three 
times that of Galactic absorption. 

4.1.2. Single DBB Fits 

A summary of fitting the single-component multicolor disc 
models from the simulations is presented in Figure QT| In 
the left-hand plot of this figure the A^h values from the DBB 
model with a fixed PO component of T=1.7 are plotted for 
disc temperatures 0.50 keV, 0.75 keV, 1 .00 ke V, 1 .50 keV and 
2.25 keV. No shading is shown for the different photon index 
slopes in this panel due to the small deviations in best-fit val- 
ues of A^h- Galactic A^h is indicated by the dashed horizontal 
line. In the right-hand panel AkT m is shown against increas- 
ing flux-ratios, where the variation of best-fit kli m for different 
values of photon index is indicated by the shaded regions. The 
T=1.7 input parameter values are shown by the solid lines. 

From these figures it can be determined that when the flux 
emission does not arise entirely from the disc component of 
the LMXB, the best-fit value of A^h from the single compo- 
nent DBB model is < Galactic absorption. In fact, even with 
a thermal flux ratio of 75% the best-fit value of absorption is 
more than 3 times smaller than that of Galactic, tending to 
as the flux ratio decreases. In the right-hand panel it is shown 
that as the input parameters of the simulations generate lower 
flux ratio spectra, so the value of IA^T; n l increases. This in- 
dicates that as the composite model becomes more dominated 
by the non-thermal component the best-fit parameters from 
the DBB model become less reliable, as was seen with the 
single-component PO model for the disc-dominant spectra. 

4.1.3. Composite PO+DBB Fits 

In addition to the single component fits, a composite PO 
and DBB model, with initial parameters set to the input val- 
ues of the simulation, was also applied to each of the spectra. 
In all cases the fit was statistically acceptable, however, from 
comparing the best-fit values to the input parameters the two 
component model was unable to adequately recover the ini- 
tial values. This is due to the degenerate nature of the disc 
and PO components of the model and indicates that, in gen- 
eral, for sources with less than 1000 counts a two component 
model will be unable to provide well constrained parameters 
(although this will be dependant on the shape of the spec- 
trum). 



4.1.4. A^h as a Discriminant 

These simulations provide a robust framework in which 
low-count single-component spectra can be used to infer 
properties of LMXBs and their spectral state. From these 
spectra it has been demonstrated that, whilst PO or DBB mod- 
els are not capable of accurately constraining the spectral pa- 
rameters of X-ray binaries when both a strong thermal and 
non-thermal component are present in the spectrum, they can 
provide important information that can be used to aid the in- 
terpretation of the source's spectral state. From the single- 
component PO model the A^h parameter is a strong indicator 
of the flux ratio between the two spectral components, where 
an elevated value of intrinsic absorption is indicative of the 
source containing a significant disc component (of at least 
25% flux for lower inner disc temperatures). In fact, this ther- 
mal to non-thermal contribution can be further constrained by 
determining the A^h absorption value from the DBB fit, where 
only sources in a thermally-dominant state provide values that 
are consistent with Galactic A^h- Further, these simulations 
have also determined the reliability of the photon index or 
ATi n values recovered from the single-component models for 
these composite sources. In the case of the PO model, when 
the source is in a hard state, indicated by a value consistent 
with Galactic A^h, the best-fit value of T provides an accu- 
rate value of the input parameter. Similarly, when the best-fit 
value of A^h from the DBB model provides values > 0, indi- 
cating the source is in a thermally dominant state, the best-fit 
ATi n provides a reliable value of the true inner disc tempera- 
ture. In all other instances neither models provide a reliable 
measure of the intrinsic source properties. 

4.1.5. Flux Determination 

In addition to the V, AT; n and A^h parameters, 0.3-8.0 keV 
values of flux from the single component models were also 
compared to the input flux values. A comparison of the de- 
rived fluxes for disc temperatures of 0.50, 1.00 and 2.25 keV, 
and photon index values 1.7 and 2.5 are shown in Figure fT2l 
In this figure the fluxes from the three different models, as 
well as a canonical PO model of Galactic Nu and T=1.7, are 
show in six panels, where the left-hand panels indicate fluxes 
for all values of ATj n from the T=1.7 models, and the right- 
hand panels presents the flux values from the T=2.5 models. 
The fluxes derived from the composite models are indicated 
by the solid black lines, the single PO model by the dot- 
dashed green lines, the best-fit values from the DBB model 
the dashed blue lines and the canonical model flux is indi- 
cated by the dotted orange line. 

From this figure it can be seen that as the disc component 
becomes more prominent in the composite model, the value of 
flux derived from the single-component PO model diverges 
from the flux determined from the composite model. The 
value of this divergence increases as the temperature of the 
inner disc decreases. Further, by comparing these changes in 
flux with Figure [9] it can be seen that in all instance where an 
elevated value of A^h is determined from the PO model, the 
intrinsic flux from the fit is significantly higher than the true 
source value, a consequence of the excess soft emission of 
the disc being modeled by the higher absorption component. 
Conversely, the intrinsic flux from the DBB best-fit model 
show a lower value than the one defined by the input model 
in all instances, up to a disc contribution of 100%. Although, 
these lower values provide a much closer value to the true 
flux than that provided by the single-component PO fit (for 
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instances where absorption column values above Galactic jVh 
are recovered). However, the flux values that are closest to the 
input model values are derived from the canonical PO model. 
This is true in most cases where the fitted A^h from the sin- 
gle PO model is significantly above the Galactic value, for 
both input T=1.7 and 2.5. However, this canonical PO model 
does provide a less accurate value of the true flux compared 
to the DBB model when a significant disc is present. There- 
fore, when providing flux values for sources with parameters 
determined from single-component models (when Nh is en- 
hanced), the flux from the canonical model should be used 
as an estimate of the true flux, except when the DBB model 
best-fit value of A^h is consistent with Galactic absorption, in- 
dicating that the source is in a TD state. In these instances 
the DBB flux value should be taken as a closer estimate of the 
source flux. 

It should be noted that the interpretation of these spectral 
models has made the assumption that these sources are not 
embedded within diffuse emission arising from the hot gas 
within the galaxy. In the instances where the galaxy contains 
a substantial ISM care should be take to determine if the en- 
hanced values of A^h are a consequence of intrinsic emission 
or caused by the surrounding medium of the source. Further, 
intrinsic absorption from the source itself, e.g. a warped ac- 
cretion disc; as proposed for the BH-GC source in NGC 4472 
by Shih et al. (2008), cannot be ruled out when enhanced 
values of A^h are determined in the spectra. 

4.2. Cool-Disc Models 

In the previous section, simulations covering the range of 
parameters that have been observed in LMXBs in three spec- 
tral states: hard, thermally dominant and the steep power law 
state, have been presented. In this section the simulations are 
extended to include spectra that have been observed in X-ray 
binaries with SPL characteristics, but have cooler inner disc 
temperatures than are typically observed for sources in this 
state (ATi n <0.4 keV). From spectral modeling of RXTE ob- 
servations of Galactic sources with this emission (e.g. XTE 
J 1550-564, Kubota & Done 2004) it has been suggested that 
these sources appear to have a cooler disc component due to 
the presence of a Comptonizing corona coupling with the disc, 
altering the intrinsic properties of the accretion disc (Done 
& Kubota 2006). This type of cool-disc spectrum has also 
been observed in ULXs, albeit with shallower photon indices 
(~1.7, e.g. Soria et al. 2007). This ULX emission has been 
likened to the spectra observed in Galactic binaries in outburst 
(Kubota & Done 2004), where the explanation of a Comp- 
tonizing component altering the properties of the accretion 
disc has again been suggested. 

To investigate the range of parameters obtained from single 
component modeling of this cool-disc spectral state, simula- 
tions were produced with input inner disc temperatures of 0. 1 , 
0.2 and 0.3 keV, and photon index values of 1 .7 and 2.5; repre- 
senting both ULX and Galactic LMXB spectra. Flux ratio val- 
ues of 75%, 60%, 50%, 40%, 25% & 0% were produced for 
each set of input parameters. A summary of the single com- 
ponent PO models of these simulated spectra are presented 
in Figure Q~3] following the same convention as in Figure [9] 
In this figure only the T=1.7 input models are indicated in 
both plots for clarity. From the left-hand panel in this figure, 
where values of Nh are presented, it can be seen that when 
the source has a high flux ratio (with fcTi„>0.1 keV) the best- 
fit value of A^h is much higher than that of Galactic, as was 
seen in the previous simulations. As the non-thermal compo- 



nent becomes more prominent the value of A^h decreases and 
in fact, for flux ratio values <25%, the best-fit value of A^h 
is either consistent with or below Galactic for disc tempera- 
tures of 0.3 keV, and below Galactic with a flux-ratio value of 
<50% for AT; n =0.2 keV. For the simuations with an inner disc 
temperature of 0. 1 keV all values of A^h are 0, only becoming 
consistent with Galactic A^h with flux ratios <25%. 

The best-fit values of A^h derived from the simulations with 
the non-thermal component described by a PO with T=2.5, 
exhibited similar variations to those presented in Figure Q~3] 
with all A^h values from the 0. 1 keV inner disc simulations 
< Galactic A^h- However, there were differences between the 
0.2 and 0.3 keV simulations from the different T values, with 
the single component fits from the steeper value of 2.5 always 
determining a best-fit A^h higher than Galactic for all flux ratio 
values, apart from 0%. 

The single-component DBB models for the cool-disc sim- 
ulations were statistically rejected in 100% of cases when the 
input photon index was 1.7. When a slope of 2.5 was used, 
the 0.2 and 0.3 keV disc temperature simulations were ade- 
quately described by single-component models in <50% of 
cases. When these fits were statistically acceptable the best- 
fit value of A^h was 0. When the inner-disc temperature was 
0.1 keV, none of the simulated spectra were statistically ac- 
ceptable. 

4.3. Different Absorption Models 

In the simulations that have been presented in this paper 
the phabs model (Balucihska-Church McCammon 1992) has 
been used to describe the absorption of X-rays in the inter- 
stellar medium. However, other models of neutral absorp- 
tion such as wabs (Morrison & McCammon 1983) and tbabs 
(Wilms et al. 2000) are also commonly used to describe the 
emission from X-ray binaries. 

To determine if the results from the phabs simulations are 
consistent with results from alternative absorption models, a 
sub-set of simulations were generated using both wabs and 
tbabs models. To investigate the range of parameters obtained 
from single component modeling of these different absorp- 
tion models, simulations of a 1000 counts were produced for 
sources with input inner disc temperatures of 0.50, 1.00 and 
2.25 keV, and with a photon index value of 1.7. Flux ratio 
values of 100%, 90%, 75%, 60%, and 40% were produced for 
each set of input parameters. A summary of the best-fit A^h 
values from the single component PO models of these simu- 
lated spectra are presented in Figure [14] following the same 
convention as in the left-hand panel of Figure [9] In this figure 
the diagonal black lines indicate the phabs model, the green 
lines the tbabs model and the blue lines the wabs model. The 
standard deviation of each derived value of A^h is indicated 
by the color-coded errorbar in the top left corner. From this 
it can be seen that all three models result in similar values 
of A^h, and are identical within the standard deviation. The 
best-fit values of A^h from the single-component DBB model 
were also compared, and all three models again resulted in 
almost identical best-fit values of A^h, with the fitted absorp- 
tion column tending to as the flux ratio decreased. These 
comparisons demonstrate that when spectrally modeling low- 
count data with single-component models, the selection of ei- 
ther wabs, phabs or tbabs in the spec tral model will result in 
the same behavior outlined in section |4~T1 

4.4. Simulations Applied to Other Galaxies 
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To investigate how the single-component best-fit parame- 
ters altered when applying the models to galaxies with dif- 
ferent values of absorption column, a sub-set of simula- 
tions were performed. For the first case, simulations were 
produced for NGC 4278, an elliptical galaxy with a line 
of sight absorption below that of the current simulations 
(1.76x 10 20 cm" 2 compared to 2.78 x 10 20 cm" 2 ) 7 . These simu- 
lations were performed using response files from spectral ex- 
traction of sources within the galaxy, and input parameters 
were set to cover inner-disc temperatures of 0.50, 1.00 and 
2.25 keV, with a fixed photon index of 1.7. Flux ratio val- 
ues from 40% to 100% were investigated, with all simulated 
spectra produced with 1000 counts to allow comparisons to 
be made with the models presented in section |4~T1 This sub- 
set of simulations was then extended to determine the best- 
fit parameters for systems with Galactic absorption values of 
5 x 10 20 cm" 2 and 1 x 10 21 cm" 2 , with input parameters defined 
to cover the same V, kT[ n and flux ratio values as for the NGC 
4278 simulations. 

A comparison of the results from these four sets of sim- 
ulations is presented in Figure Q3] where the excess absorp- 
tion column values from the single-component PO models 
are shown. These excess A^h values are defined as the best-fit 
value of A^h normalized by the simulated Galactic absorption 
(A^ho)- In this figure the black lines indicate the simulations 
produced for NGC 3379, the orange lines the excess A^h val- 
ues of the NGC 4278 simulations, blue the Af H o=5 x 10 20 cm" 2 
simulations and the green lines the Nno=l x 10 2I cm~ 2 simula- 
tions. In the top left corner color-coded error bars indicate the 
standard deviation of each set of simulations. What can be 
seen from this figure is that all simulations exhibit the same 
behavior of increasing values of excess A^h as the thermal 
component of the source becomes more prominent. In fact, all 
excess absorption values from the four different sets of simu- 
lations are identical within the standard deviations. The spec- 
tra from these simulations were also fitted to single compo- 
nent DBB models, where the best-fit A^h values followed the 
behavior shown in Figure QT| .This indicates that these simu- 
lations are applicable in all galaxies, where single-component 
best-fit values can be used to determine the presence of ab- 
sence of a thermal component, as well as the intrinsic source 
parameters. 

4.5. Decision Tree 

Figure [16] provides a summary indicating how the best-fit 
parameters derived from the single-component models can be 
used to determine the spectral state of a source. This is pre- 
sented in the form of a decision tree, indicating how best-fit 
values of A^h from the PO and DBB models can be used to dis- 
criminate between NT, TD and intermediate states, including 
sources containing a cool-disc component. 

5. DISCUSSION 

From the spectral analysis of the eight sources presented 
in this paper a diverse range of temporal characteristics have 
been identified in the bright (L x > 1.2 x 10 38 erg s" 1 ) LMXB 
population of NGC 3379. Here the properties of the four 
sources determined to exhibit spectral variability are pre- 
sented individually and their best-fit parameters are compared 
to the spectral simulations discussed in section [4] The spec- 
tral parameters of the remaining four sources are also pre- 

7 This galaxy is the subject of a companion paper Fabbiano et al. (2010). 



sented and the overall properties of this bright population are 
discussed. 

5.1. A Black Hole GC-LMXB 

Of the three GC-LMXBs that are presented in this paper, 
S42 is both the most luminous and the most variable. From 
Figure |2]it can be seen that there is a significant increase (and 
subsequent decrease) in flux between observations implying 
AL X ~ 2.5 x 10 38 erg s -1 in the simple assumption of a canoni- 
cal PO spectrum used to derive Figure|2](B08). This pattern of 
variability is fairly common for LMXBs (van der Klis 1994). 
The large flux variation observed very strongly indicates that 
S42 is not a superposition of multiple NS-LMXBs but a sin- 
gle source. In fact, if it is conservatively assumed that the 
total luminosity is from the contribution of multiple sources, 
the luminosity from a single source would have a lower limit 
of ~ 2.5 x 10 38 erg s" 1 (increasing to ~ 7 x 10 38 erg s" 1 if 
the emission arises from a single source). This luminosity is 
above the Eddington limit for a neutron star binary (although 
would be close to the limit of a heavy (2-3 M Q ; Kalogera & 
Baym 1996) NS, or a 1.4 M NS with a He or C/O donor) 
and indicates that the source is likely to be a black hole binary 
system. 

Based on the X-ray luminosity, HR values and colors of this 
source, the spectra were separated into three different groups. 
In all three cases both single component models provided sta- 
tistically acceptable fits to the data. In the PO model the 
spectral groups provided a range in T of ~ 2.1 - 1.7, which 
is within the range of a hard state BH-LMXB. However, in all 
three cases a large intrinsic Nn, of order 6 times the Galac- 
tic column, was required to provide an adequate description 
of the data. Such large Nu columns are not usually seen 
in LMXBs and would be especially difficult to explain in a 
sub-Eddington low-state source, where outflows are not ex- 
pected. One explanation could be that the lower luminos- 
ity state is a consequence of foreground absorption from a 
warped disc (as advocated for the variable BH-GC source in 
NGC 4472, which showed variability consistent with arising 
from changes in Nu; Shih et al. 2008). However, since there 
is little variability detected in Nu between spectral groups for 
S42, this mechanism is unlikely to be the cause of the change 
in Lx- Instead, from the simulations presented in section [4] 
there is an indication that when a source has a significant 
flux contribution from a thermal disc component, the single- 
component PO model will provide an elevated absorption col- 
umn, caused by the soft emission arising from the disc. Fur- 
ther, the DBB single-component modeling from the simula- 
tions indicates that when a source is in a TD state the best-fit 
absorption column >0, becoming consistent with Galactic ab- 
sorption as the accretion disc completely dominates the flux 
contribution (see Figure fTTTl. 

Indeed, the DBB best-fit models of S42 suggests that this 
source is predominately in a TD state over observations 1-4, 
where Nu ranges between 0.6-3.2 x 10 20 cm" 2 , possibly enter- 
ing an intermediate state in observation 5, where its likely 
that the thermal and non-thermal flux contribute equally to the 
emission. This disc component in the TD state appears to be 
at its most dominant when the source is at its most luminous, 
in spectral grouping 2, where the best-fit value of A'h is con- 
sistent with Galactic and kT m ^l.5 keV, which infers a mass 
of -4 M (cf. Gierlihski & Done 2004). This best-fit model 
provides a luminosity of ~ 9.1 x 10 38 erg s , which further 
strengthens the interpretation that this source is a black-hole 
binary, where flux variability indicates that there must be a 
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source with Lx> 4 x 10 38 erg s -1 within the GC. The increase 
in temperature coupled with the increase in luminosity ob- 
served between the spectral groupings could be explained as 
the disc's response to a temporary increase of the accretion 
rate, or could be due to the failure of the single-component 
model to adequately describe the properties of the spectra as 
the PO component becomes more prominent. However, the 
statistics in these spectra do not constrain more complex mod- 
els (e.g. PO+DBB) to allow this to be investigated further. 

What is very likely from these spectral properties is that a 
single black-hole binary has been detected within a globular 
cluster in NGC 3379, in addition to the BH-GC within NGC 
4472 (Maccarone et al. 2007) and the suggested IMBH-GC 
in NGC 1399 (Irwin et al. 2010). Further, it has also been 
determined that this source is in a thermally dominant state 
when it is at its most luminous. Maccarone et al. (2007) find 
that the BH-GC in NGC 4472 is in a very luminous (~7x 
10 5 Lq) and metal poor GC and Irwin et al. (2010) determine 
that the source in NGC 1399 is in a luminous (~3x 1O 5 L ) 
but metal rich GC. In NGC 3379 not only S42, but also the 
two other GC-LMXB sources; S41 and S67, are also found in 
luminous GCs (B08), with luminosities of order 1.5x 10 5 L Q , 
3x 1O 5 L and 2x 10 5 L Q respectively, therefore following the 
trend identified in previous GC-LMXB studies (e.g. Kundu 
et al. 2007). Using the color-metallicity relations in Smits et 
al. (2006), and the GC colors reported in B08, we find that 
S41 and S67 are in 'blue' low-metallicity clusters (V-I=0.85 
and V-I=0.87; ~ 1/10 solar), whereas S42 is in a 'red' high 
metallicity cluster (V— 1=1.12). 

These comparisons, albeit based on only five objects, are in 
agreement with the conclusion that GC mass may be a factor 
in retaining a BH-LMXB. However, given that this paper only 
concentrates on the brightest LMXBs within NGC 3379, we 
do not speculate on the suggestion that the more massive GCs 
host the more luminous LMXBs (Kundu et al. 2007). In 
fact, from B08 it has been determined that there are two less 
luminous LMXBs that reside in more luminous GCs than the 
sources presented here. Further, in the cases of S41, S67 and 
the source in NGC 1399, without large temporal luminosity 
variations it cannot be ruled out that these GCs may contain a 
superposition of less luminous NS rather than a single bright 
BH. However, due to the large variability in both S42 and the 
source in NGC 4472, it can be confirmed that these sources 
are bright single objects. 

5.2. A Field BH-LMXBs 

S86, the only confirmed field source with Lx> 10 38 erg s -1 
within NGC 3379, exhibits similar Lx variability to that of 
the BH-GC, S42. However, in this case, when the source is 
observed in its less luminous state the spectra are consistent 
with either a hard state, where T ~1.5 and no extra absorp- 
tion is required, or possibly a cool-disc model. This latter 
interpretation is proposed as the best-fit value of A^h tends 
to when left free to vary and the single-component DBB 
model is statistically rejected when applied to these spectra. 
When investigating the cool-disc simulations, such behavior 
was observed in both these single-component model fits (see 
section l4!2l . In the higher flux observations the DBB model 
describes the spectra well, where the best-fit inner disc tem- 
perature is ~ 1 .4 ke V, consistent with the TD state of a ~5 
M BH. The PO model also provides an acceptable descrip- 
tion of the spectra in this higher flux state, with r ~1.6, typ- 
ical of hard state spectra. However, in this more luminous 
state, whilst the spectral index indicates a NT source, the ab- 



sorption column is a factor of 4 greater than the Galactic A^h- 
This could again point to the single component model failing 
to describe the nature of the accretion, indicating a change in 
spectral state from the NT to the TD state (as investigated in 
the simulations in section or it could indicate the presence 
of a warped disc, as discussed in Shih et al. (2008). In both 
models the source luminosity exceeds the NS Eddington lu- 
minosity, with Lx ranging between ~ 7 x 10 38 erg s" 1 in the 
NT state (or cool-disc) observations (where this high Lx also 
indicates that this source is not residing in the NT state), to 
~ 5 x 10 38 erg s" 1 in the TD state spectra, which is therefore 
consistent with emission from a stellar-mass BH. 

5.3. Spectral State Change of Central Region LMXB 

S74 resides in the central region of the galaxy where the 
optical data become confused and therefore has an unclassi- 
fied optical correlation. The spectral behavior of this source is 
very stable in the first four observations, where the spectra are 
well fitted with a T ~1.6 PO model with no intrinsic absorp- 
tion, consistent with a NS or BH-LMXB in the hard-state. Al- 
though, this source could have a cool-disc component, given 
that A^h tends to when left free to vary in the PO model, 
and the DBB model provides a poor description of the data. 
In observation 5 the spectral emission changes dramatically, 
increasing in flux and altering in spectral state, where the PO 
model indicates that the source has become softer and the ab- 
sorption column has increased. This is again interpreted as 
an indication that the disc emission has become more promi- 
nent in the source, where the DBB model describes the spec- 
tra well, although A^h left free to vary provides a best-fit value 
well below Galactic and has therefore been frozen in the fi- 
nal model. These parameters indicate that the source is in 
a TD state, with an inner disc temperature of ^0.7 keV and 
an intrinsic luminosity of ~ 7.5 x 10 38 erg s" 1 , that, if follow- 
ing the M oc T" 4 relation, indicates that the source is a ^75 
M© black hole. However, as the absorption column from the 
DBB model is below the Galactic value it is possible that a 
non-thermal component is present in the source, therefore this 
derived mass should be treated as an upperlimit. 

5.4. A Luminous Cool-Disc Source 

The final source with spectral variability is S102, which 
is just outside the HST FOV. From the plots indicating the 
variability in Lx and HR (Figure |2} alone it appears that this 
source is gradually transitioning from a bright hard source to 
a lower luminosity soft source. However, from the spectral 
groupings it has been determined that the source actually re- 
sides in distinct states. In the first observation the spectrum is 
well described by a steep power law {V ^2.6) with no intrin- 
sic absorption. This photon index value is steeper than would 
be expected for a source in the hard-state (RM06: T <2.1), 
suggesting that the source is in a SPL state. However, given 
that the single-component DBB model is statistically rejected 
when fitted to this spectra, this source could be in an interme- 
diate state with a cool-disc component (see section l4~2l i. Al- 
though, within errors, both A^h and V are consistent with the 
expected values of a SPL state source, and there are too few 
counts in the data to allow additional components to meaning- 
fully constrain the spectral parameters. 

In observations 2 & 3 the spectra exhibit very similar prop- 
erties to those observed in observation 1 although, with a 
greater number of counts, it is clear that the fit value of A^h 
from the PO model is significantly below that of Galactic ab- 
sorption (see Figure [7} and consequently a two-component 
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model has been applied to the spectra. From this best-fit 
model the thermal emission has been determined to have an 
inner-disc temperature of fcT; n — 0.14 keV, with a photon in- 
dex of 1.6. A large absorption column is required for this fit, 
although this is not well constrained and is consistent with 
Galactic absorption (see contour plots in Figure [8j. Such a 
spectrum is consistent with the properties observed in some 
ULXs (e.g Soria et al. 2007) and indeed, the intrinsic lumi- 
nosity of SI 02 in observations 2 & 3 is — l.lxl0 39 erg s" 1 , 
which does fall under the classification of a ULX, albeit at 
the very low end of the luminosity range. Further, the flux ra- 
tio from both of these observations indicates that 40% of the 
emission arises from the disc-component, slightly higher than 
the contribution determined in Soria et al. (2007), where the 
thermal emission contributed 25% of the observed flux. From 
the M ocT~ 4 relation the temperature of this fit implies an in- 
termediate mass black hole of ~ 5 x 10 4 M Q , although, as dis- 
cussed in section l4~2l this disc may appears cooler due to the 
prescence of a Comptonizing corona altering the properties of 
the disc (Done & Kubota 2006) and the true mass of the black 
hole may be much lower. 

In the final two observations of this source the luminos- 
ity decreases significantly and the flux is well described by a 
single-component PO model, with A^h consistent with Galac- 
tic absorption and T — 1.85, indicating that the source is likely 
to be emitting in a NT state. This source therefore provides 
observational information of a black-hole binary source tran- 
sitioning from a steep power law cool-disc state to a hard 
state, where the luminosity falls from — l.lxl0 39 erg s _1 to 
-3.8 x lO^ergs" 1 . 

5.5. The Overall Population 

In addition to the four sources that exhibit spectral variabil- 
ity, the four other sources included in this paper have also 
been analysed. From their spectral parameters, two are con- 
sistent with being in an NT state (S41 & S77) with photon 
index value of —1.9 and —2.0. The two other sources (S67 & 
S103) have large values of intrinsic A^h from the PO model, 
indicating that a significant disc is present in the spectra. In 
the case of S67, the DBB model provides a value of A^h con- 
sistent with Galactic absorption, suggesting that the source is 
in a TD state, whilst SI 03 has a value tending to and is 
therefore likely to be in an intermediate state with flux con- 
tributions from both the thermal and non-thermal components 
of the binary. 

In Figure [T7] the spectral properties of the binaries investi- 
gated in this paper are presented. In this figure the spectral 
groupings of each source have either been plotted in the left 
or right panel, depending on their inferred state, where in the 
left panel log kT\ n against log Lx is presented and the right 
panel shows T against log Lx- When sources are likely to be 
in a spectral state that is not TD (or NT) the errors bars are 
presented as dotted lines. The best-fit values from the PO and 
DBB components of the luminous, cool-disc observations of 
S102 are indicated by the dashed error bars. In the case of the 
TD state for this source, this point is presented as an insert 
within the main plot, where a lower disc temperature range 
is presented. In addition to the sources, the L-T relations 
for a non-rotating BH (cf. Gierlihski & Done 2004; eq. 3, 
assuming an inclination angle of 60 degrees and / co i=l-8) of 
masses 5M Q , 1OM and 2OM are indicated by the dashed 
lines in the left-hand plot, and an additional line indicating 
the L-T relation for a BH mass of 500M Q is included in the 
insert plot. In both panels, sources are indicated by symbol 



and color where labels are given in the bottom right corner. 

From comparing both panels it can be seen that, in general, 
when a source is emitting at a lower luminosity, the spectrum 
indicates that the source is in a NT state, dominated by the PO 
component. This behavior is commonly seen in Galactic bi- 
naries, where this state was previously defined as the low-hard 
state, denoting both the photon index and flux of the source. 
However, in McClintock & Remillard (2006) it was advocated 
that luminosity should not be part of the classification scheme, 
as, while many Galactic binaries do become harder as their lu- 
minosity decreases, there are counterexamples indicating that 
this behaviour is not ubiquitous (e.g. GX 339-4; RM06). In- 
deed, even within the sample of sources presented here there 
is an indication that we observe a source transitioning between 
a TD and NT state whilst exhibiting very little change in lumi- 
nosity (S86), although it is possible that these NT states (obs 
1, 3 & 5) are actually in a cool-disc state. The other bright 
sources in the Lx - T plot are observations 1 and observations 
2&3 of S102. From spectral modeling it has been determined 
that this source is actually in a steep power law or interme- 
diate cool-disc state in observation 1 and in an intermediate 
cool-disc state in observations 2&3, both of which typically 
result in a higher level of flux than the NT state. 

In the left-hand panel the five sources determined to be in a 
TD or intermediate/TD state are presented (along with an in- 
sert presenting the TD properties of SI 02 in a luminous cool- 
disc state). The five sources within the main plot all reside in 
the same L-T parameter space as Galactic black-hole bina- 
ries (cf. Miller et al. 2004; Kajava & Poutanen 2009), albeit 
at the brighter end of the relation. S102, on the other hand, 
exhibits spectral properties that do not follow the typical rela- 
tions observed in Galactic binaries (e.g. Miller et al. 2003; 
Stobbart et al. 2006) and could sugg est the presence of an 
IMBH. However, as discussed in {35.41 with such a cool-disc, 
while it is consistent with this interpretation, it is not proof of 
an IMBH. 

In this figure S42 is the only LMXB that has been deter- 
mined to exhibit multiple spectral properties whilst in the TD 
state and is denoted here by the open squares. By comparing 
these points with the dashed lines on this plot, which indicate 
the LocT 4 relation; the expected signature for optically-thick, 
geometrically thin accretion discs (Done, Gierlihski, Kubota 
2007), it appears that this source does not follow the stan- 
dard relation, but instead roug hly follows Lcxr 125 . A flat- 
ter relation of LocT 2 has been observed in other black hole 
binary systems at higher temperatures (—1 keV) (Kubota & 
Makishima 2004) which has been explained as the transition 
of the disc from a standard accretion disc (Shakura & Sun- 
yaev 1973) to a slim disc solution (Abramowicz et al. 1988) 
as the accretion rate increases. However, Gierlihski & Done 
(2004) argue that such a transition should only occur above 
^Edd, which is greater than the luminosity of sources — 1 keV 
in Kubota & Makishima (2004), and they instead suggest that 
this flattening is associated with a change in the color temper- 
ature correction (/ co i)- 

In the case of S42, the inferred mass of the BH from the 
M-T relation is — 4M , therefore the observed luminosity 
of 9.1xl0 38 erg s" 1 would require super-Eddington accretion 
for a BH of this mass. This suggests that the flatter relation 
in Figure [17] is a consequence of the source transitioning to a 
slim disc, although a change in the color temperature correc- 
tion cannot be ruled out. However, this suggestion of a flatter 
relation should be treated with caution due to, not only the 
limited data points, one of which is from the source in an in- 
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termediate state (observation 5) and therefore is likely to have 
an incorrect kT m , but also due to the uncertainty in Lx from 
the single-component spectral modeling. 

6. CONCLUSIONS 

In this paper we have investigated the spectral properties of 
eight bright (L x > 1.2x 38 erg s" 1 ) sources within NGC 3379, 
selected from the B08 catalog. From deep multi-pointing 
Chandra observations we have been able to determine the 
source properties at different epochs, and thereby characterize 
source spectral variability. A set of simulations used to infer 
the spectral states of these sources has also been presented. 

Due to the luminosity cut imposed on this sample, all of the 
eight sources are either NS emitting at, or above the Edding- 
ton luminosity, or BHBs, with source luminosities ranging be- 
tween 1 .6- 10.8 x 10 38 erg s" 1 . From optical correlations three 
sources have been determined to be BH-GCs, one is detected 
in the field, while the remaining four sources have insufficient 
optical data to classify a correlation. In all cases (bar S102) 
single component PO and DBB models have been applied to 
the data and have provided statistically acceptable fits. 

To aid the interpretation of these single component mod- 
els, simulations of low-count (250, 500 and 1000) data have 
been presented. This work demonstrates that it is possible to 
determine if a source is in a hard, thermally dominant, or in- 
termediate state by comparing the best-fit values of A^h from 
the two spectral models. In cases where the line-of-sight ab- 
sorption in the PO fit is significantly higher than that of the 
Galactic value it is likely that a strong thermal component is 
present. If it is then determined from the DBB model, that 
this same source has a best-fit A^h consistent with the Galactic 
value it points to the source being in a TD state. If instead the 
DBB best-fit A^h tends to 0, it is likely that the source has both 
strong thermal and non-thermal components and can be said 
to be in an intermediate state. Although we add the caveat that 
whilst an elevated value of A^h from the single-component PO 
model is indicative of a thermally dominant state, intrinsic ab- 
sorption, arising from the disc could also explain this result in 
some cases. We have also provided simulations of cool-disc 
spectra, allowing us to identify when a source is in a high 
accretion state with a cool-disc component. A summary of 
how the best-fit values from single-component models can be 
interpreted is provided in Figure [T6l 

From these simulations it has been determined that when 
a source has prominent disc emission as well as a signif- 
icant non-thermal component, the luminosity derived from 
the single-component DBB model should be taken as a lower 
limit of Lx- The value inferred from the single-component 
PO model will be much greater than the true source value 
and should be disregarded. However, when a canonical PO 
model (r=1.7 and A^h frozen to the Galactic absorption value) 
is applied to the data, the derived Lx value will provide the 
closest value to the true source luminosity. When the source 
is in a TD state the X-ray luminosity from the DBB should be 
adopted. 

We have also performed these simulation for a range of dif- 
ferent absorption values, confirming that the results are not 
significantly dependant on the value of Galactic absorption, 
demonstrating that these simulations can be used to interpret 
a variety of extra-galactic binary spectra. Further, different 
absorption models have also been been used in a sub-set of 
simulations, and these results have confirmed that with low- 
count data the choice of absorption model does not affect the 
interpretation of spectral state. 



Four of the eight sources presented in this paper have been 
determined to exhibit spectral variability and observations 
have been subsequently been grouped to reflect these varia- 
tions. By comparing the spectral properties of these group- 
ings to our simulations, spectral transitions have been iden- 
tified. The most luminous source detected within a GC, 
S42, has been determined to show significant variability (Lx> 
4 x 10 38 erg s" 1 ) and can therefore be confirmed as a strong 
BHB candidate, predominantly residing in the thermally dom- 
inant state. The best-fit inner-disc temperature of this source 
is 1.5 keV, and we estimate a black hole mass of ^4M , pro- 
viding further evidence (in addition to Maccarone et al. 2007 
and Irwin et al. 2010) that GCs do retain black hole bina- 
ries. If we assume that the two further GC-LMXB sources in 
this paper are also individual black holes, this indicates that 
30% of the GC-LMXB correlations in NGC 3379 contain a 
black hole binary, which constitutes a significant fraction of 
the GC-LMXB population. 

A transition from a luminous steep power law state to a 
hard state has been identified in SI 02 where we have iden- 
tified a binary with a cool-disc component (fcTj n =0.14 keV 
and r=1.6) emitting around the Eddington luminosity of a 
1OM black hole (Lx^ 1 x 10 39 erg s" 1 ). Taken at face value, 
the temperature of this cool-disc could suggest an IMBH. 
This source then declines over the three months between 
pointings to a source in a hard state, with T ~ 1.85 and 
Lx=3.8 x 10 38 erg s -1 , consistent with NT state Galactic black 
hole binaries (RM06). The bright state of this source is sim- 
ilar to spectra observed in some ULXs (Soria et al. 2007), 
and indicates the similarities between these bright sources and 
more 'normal' stellar mass black hole binaries in high accre- 
tion states. 

The spectral properties of the eight sources are largely con- 
sistent with the parameters that have been observed in Galac- 
tic LMXBs, with sources in a NT state exhibiting a range of 
r=1.5-1.9 and sources in a TD state with /tT in =0.7-1.55 keV. 
We have also shown that this population is consistent with 
the general trend of increasing luminosity as sources become 
softer, transitioning from a NT to TD state. The L—T relation 
has been investigated for all the sources in a TD state, with 
the parameter space of these LMXBs (excepting the luminous 
cool-disc state of SI 02) being consistent with Galactic obser- 
vations, albeit at the bright end of the L oc T 4 relation. S42, 
the BH-GC, was the only source to be observed in the TD 
state in multiple pointings. The L—T relation of this object 
is flatter than the typical relation, with luminosity increasing 
with temperature as ~ T L25 . Flattening has been observed in 
Galactic black hole binaries (Kubota & Makishima 2004) and 
it has been suggested that this could be due to the standard 
Shakura-Sunyaev accretion disc evolving into a slim disc, al- 
though this flattening could also be explained by a change in 
the disc color temperature relation. 
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TABLE 1 

Observation Log and Source Counts 





Obs 1 


Obs 2 


Obs 3 


Obs 4 


Obs 5 


All 


Opt. corr 


OB SID 


1587 


7073 


7074 


7075 


7076 






Exposure (ksec) 


29.0 


80.3 


66.7 


79.6 


68.7 


324.2 




Net Counts 


S41 


7.5±10.0 


146.0±13.5 


143.6±13.2 


150.0±13.6 


200.0±15.4 


718.7 ±28.4 


GC 


S42 


170.4±14.1 


457.8±22.5 


269.0±17.6 


512.8±23.8 


331.5±19.4 


1741.8±43.0 


GC 


S67 


95 .4± 11.1 


239.5±16.9 


219.0±16.1 


243.2±16.9 


198.7±20.9 


1008.9±33.4 


GC 


S74 


55.8±9.0 


181.1±14.8 


147.5±13.7 


158.2±14.0 


580.7±25.5 


1165.4±35.9 




S77 


25.5±6.7 


111.1±12.2 


90.0±11.2 


92.0±11.3 


71.7±10.1 


414.4±22.5 




S86 


195.7±15.2 


572.9±25.3 


371.4±20.6 


620.4±26.2 


352.4±20.2 


2114.3±47.7 


F 


SI 02 


165.5±14.0 


459.4±22.6 


332.5±19.5 


288.3±18.2 


140.2±13.2 


1389.0±38.7 




S103 


75.1±10.6 


131.3±12.6 


127.1±13.4 


170.7±15.1 


143.1±14.1 


636.9±26.6 





NOTE. — Optical correlations: GC indicates an association with a globular cluster, F indicates a field source, which is within the HST FOV and has no optical association. - 
indicates that there is insufficient optical coverage for this source to determine an optical association. All of these classification are fully described in B08, §2.6. 



TABLE 2 

Summary of Long- and Short-term Source Variability 



Source No. Short-term variability Long-term variability 





Obs 1 


Obs 2 


Obs 3 


Obs 4 Obs 5 Var. 


Sign. 


S41 








V 


3.9cr 


S42 








V 


5.9cr 


S67 








N 


1.0(7 


S74 


P 






V 


15.8c- 


S77 








V 


2.4cr 


S86 




P 




V 


6.0a 


S102 








V 


10.7cr 


SI 03 






P 


V 


2.0cr 



NOTE. — Individual variability is based on the Kolmogorov-Smirnov test: - indicates no variability, V, indicates variable source (with values >99% confidence) and P indicates 
possible variable sources (with variability values >90% confidence). The long-term variability from all five observations was defined by the chi-squared test where V indicates 
variability, and N indicates a non-variable source. Sign indicates the significance in change in luminosity between the highest and lowest flux observations. All of these classification 
are fully described in B08, §2.4. 
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TABLE 3 

Summary of the best-fit parameters of the source spectra. Errors are 

GIVEN TO lcr 



Source 


Obs. No. 


Net Counts 


Model 


X> 


P x 2 


N H 
xl0 iu 


Parameters 

r 




X 1() 38 erg s' 1 
(0.3-8.0 keV) 


L x Range 


S41 GC 


All 


719 


PO 


25.89/24 


0.36 


2.8Fr 


i m +012 
iv ^-0.11 


- 


2.87 


1.99-4.51 


S42 GC 


1&3 
1&3 


439 
439 


PO 
DBB 


10.32/15 
12.00/16 


0.79 
0.74 


I7.8£g 

2.8Fr 


2.09«;J0 


0.89^ 


5.60 
4.78 


4.21-8.15 
3.18-6.08 




1&3 


439 


DBB 


11.32/15 


0.73 




- 


4.15 


3.58-6.55 




2&4 
2&4 


971 
971 


PO 
DBB 


32.86/40 
40.05/40 


0.78 
0.47 


18.8!" 

3.2!|;j 




- 

1 --' u -0.13 


12.30 
9.10 


12.22-12.38 
7.79-10.76 




5 
5 


332 
332 


PO 
DBB 


3.85/12 
5.44/13 


0.98 
0.96 


12 l+ 4 - 2 
1Z.1_ 7 5 

2.8Fr 


1 75+ 013 
1 ■ ' J -0. 13 


1 - zl -0.17 


8.26 
6.14 


7.80-9.60 
3.99-6.94 




All 
All 


1002 
1002 


PO 
DBB 


31.50/34 
29.92/35 


0.59 
0.71 


17 4+7.3 
2.8Fr 


1 57 +0 - 17 


1 55 +<m 
55+8:18 

1 - JJ -0.17 


4.48 
4.19 


2.96-5.65 
2.71-5.16 




All 


1002 


DBB 


29.83/34 


0.67 






4.67 


9 S2-^ 09 


o /4 


1 91 jPr/1 




rU 


1 A 99/1 A 
10.ZZ710 


U.4J 


z.orr 


1 fit +0.15 




9 "7/1 

Z. /4 


l.J 1— J.JJ 




123&4 


553 


DBB 


26.39/16 


0.05 


2.8Fr 


1 01 +012 


1.95 


0.97-2.38 




C 
J 

5 


555 


rU 

DBB 


JJ.JJ/ZL 

35.68/22 


U.UJ 
0.03 


1 i £+3.5 
2.8Fr 


9 iq+0.11 
2 - 13 -0.14 




iz. jy 
7.65 


1 1 1 Q "7/1 

6.63-7.89 


S77 


All 


368 


PO 
DBB 


6.72/7 
7.61/8 


0.46 
0.47 


2.8Fr 


2 03+ ' 25 


_ 

0-98^3 


1.60 
1.28 


0.37-2.20 
0.49-1.71 


S86 F 


13&5 


920 


PO 


35.50/37 


0.54 


2.8Fr 


1 50+ 009 
■ J -0.09 


- 


7.07 


6.15-7.95 




2&4 
2&4 


1193 
1193 


PO 
DBB 


39 38/54 
42.17/51 


0.86 
0.81 


13 1 +2 ' 4 
2.8Fr 


1 65 +01 ° 

luJ -0.08 


1 35+ - 07 
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NOTE. — In column 1 optical correlations (or lack thereof) are indicated by the symbols: GC indicating a GC-LMXB source and F indicating that the LMXB has been confirmed 
as a field source. No symbol indicates that there are insufficient data to establish the presence of an optical counterpart. When the Cash statistic has been used in preference to x 2 
a C is included in columns 5 & 6. A value of 2.8Fr in column 7 denotes that the best-fit model has had the absorption column frozen to the Galactic value. Column 1 1 presents the 
range in X-ray luminosity where the lowest value indicates the la lower-bound value from the lowest flux individual observation from the joint fit. The higher value indicates the la 
upper-bound value for the most luminous individual observation from the joint fit. For best-fits from individual observations the range indicates the lower- and upper-bound values 
from that single observation. 
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FIG. 1 . — Adaptively smoothed full-band (0.3-8.0 keV) X-ray image of NGC 3379. Overlaid in this figure are labeled circular regions indicating the eight 
source presented in this paper and HST FOV. 
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FIG. 2. — Plots of the eight bright sources in NGC 3379, adapted from B08. Each plot summarizes the variations in properties of each source between each 
observation. In the main panel the long-term light curves are shown. In the second panel down, the hardness ratios are indicated. These are defined to be HR 
= H— S / H+S, where H is the number of counts in the hard band (2.0—8.0 keV) and S is the number of counts in the soft band (0.5-2.0 keV). In the third and 
fourth panels the color ratios; C21 and C32, are plotted, where C21=/ogS2+/ogSl and C32=-/ogH+/ogS2. For the color ratios the bandwidths are defined to be 
S 1=0.3-0.9 keV, S2=0.9-2.5 keV and H=2.5-8.0 keV. From B08 it was determined that all sources apart from S67 display long-term flux variability. 
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FIG. 3. — Contours for PO (left) and DBB fits (right) for S42. In both figures the black solid line represents the 2<r contour for observations 1 & 3, the red 
dashed line observations 2 & 4 and the blue dotted line observation 5. Best-fit values are indicated by the cross hair within each contour (these values differ from 
the best-fit parameters presented in Table[3]when Nu is below Galactic absorption). In some instances cross hairs are not readily observable, in these cases the 
best-fit value of Nu tended to when left free to vary. The solid vertical line indicates Galactic Nr. From the simulations presented in section [4] the spectral 
properties of S42 suggest that this source is predominately in a thermally dominant state, entering an intermediate state in observation 5 (see text for details). 
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FIG. 4. — Short-term lightcurve of S86 from observation 2. The x-axis provides time in seconds, beginning from the start of the observation, and the y-axis 
provides count-rate in count s between 0.3—8.0 keV. Binning of 10000 sec was used when extracting counts to provide adequately constrained count-rates to 
identify variability. 
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FIG. 5. — Contours for PO (left) and DBB fits (right) for S86. In both figures the black solid line represents the 2cr contour for observations 1,3 & 5, the red 
dotted line observations 2 & 4. Best-fit values are indicated by the cross hair within each contour (these values differ from the best-fit parameters presented in 
Table[3]when Nu is below Galactic absorption). In some instances cross hairs are not readily observable, in these cases the best-fit value of Nu tended to when 
left free to vary. The solid vertical line indicates Galactic Nu. From the simulations presented in section|4]the spectral properties of S86 indicate that this source 
is transitioning between a hard state (or possibly a cool-disc state) to a thermally dominant state (see text for details). 
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FIG. 6. — Contours for PO (left) and DBB fits (right) for S74. In both figures the black solid line represents the 2cr contour for observations 1,2,3 & 4, the red 
dotted line observations 5. Best-fit values are indicated by the cross hair within each contour (these values differ from the best-fit parameters presented in Table 
[3]when Nn is below Galactic absorption). In some instances cross hairs are not readily observable, in these cases the best-fit value of N\{ tended to when left 
free to vary. The solid vertical line indicates Galactic N^. Based on the simulations presented in section|4] S74 is likely to be in a hard state in observations 1—4, 
transitioning to a thermally dominant state in observation 5 (see text for details). 
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FIG. 7. — lcr and 2tr contours for the PO best-fit for S102. The solid black lines represents contours for observation 1, the red dashed lines show contours 
for observations 2 & 3 and the dotted blue lines observations 4 & 5. Best-fit values are indicated by the cross hair within each contour (these values differ from 
the best-fit parameters presented in Table[3]when Nu is below Galactic absorption). In some instances cross hairs are not readily observable, in these cases the 
best-fit value of N\\ tended to when left free to vary. The solid vertical line indicates Galactic N\\. From the simulations presented in section [4] SI 02 is likely 
to initially be in a steep power law state, transitioning to a luminous cool-disc state in observations 2 & 3 and then changing to a hard state in observations 4 & 5 
(see text for details). 
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FIG. 8. — lcr and 2a contours for the two-component PO+DBB modeling of observations 2&3 of S102. Left indicates values of Nn and jfcTj„ and the right 
panel presents the T and kT m values from the same fit. The best-fit values are indicated by the cross hairs in both panels. In the left panel the solid vertical line 
indicates Galactic N\\ . 
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Flux Ratio Flux Ratio 



FIG. 9. — Results of the best-fit values derived from the single-component PO model applied to the 1000 count simulated spectra (detailed in section[4}, where 
all spectra were simulated with Galactic Nu. The left panel indicates the values of Nu plotted against the flux ratio of the disc-to-total flux. The right panel shows 
the difference (A) between the input and best-fit value of T. In both plots the solid diagonal lines label input V = 1.7 for each simulation, for input &T; n values of 
0.50, 0.75, 1.00, 1.50, and 2.25 keV. The shading indicates the spread in best-fit values for the full range of photon indices 1.5-2.5. In the left panel the horizontal 
dashed line indicates the Galactic Nu an and the vertical solid lines are grid line. The error bars in the upper-left comer indicate the standard deviation of the 
0.75 (large) and 0.25 (small) value of Nu. In the right panel the kT m values are presented in different colors for clarity. It should be noted that different shading 
patterns of green, blue and gray are a consequence of overlapping shading and do not indicate additional information within the figure. 




0.2 0.4 0.6 0.8 1 

Flux Ratio 



FIG. 10. — This figure presents the same information as is shown in the left panel of Figure |9]but also includes the best-fit values of Nu derived from lower 
count simulations of input parameters T=1.7 & 2.5 and &Tj n =1.00. The 500 count data is plotted in blue and the 250 count data in green. The results of the 
r=1.7 simulations are labeled by the solid diagonal lines and the best-fit values of Nu from the T=2.5 simulations for 60%, 75%, 90%, and 100% disc-to-total 
flux ratios, are indicated by the outer edge of the shading (follwing the findings of the 1000 counts simulations, intermediate values of V are expected to reside 
within the shaded region). Error bars in the top left indicate the standard deviation for 0.75 (large) and 0.25 (small) flux ratios, these are also color-coded. The 
horizontal dashed line indicates Galactic Nu . 
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FIG. 1 1 . — This figure presents the best-fit values derived from the single-component DBB model applied to the simulated spectra with input values of Nn set 
to Galactic absorption (detailed in section[4}. As in Figure|9] the left panel indicates the values of Nu plotted against the flux ratio of the disc-to-total flux and the 
right panel shows the difference (A) between the input and best-fit value of kl m . In both plots the solid diagonal lines labels input T = 1.7 for each simulation, 
for input KT m values of 0.50, 0.75, 1.00, 1.50, and 2.25 keV. Shading is not presented in the left panel due to the small variations of output Nn values but in the 
right panel it indicates the spread in best-fit values for the full range of photon indices 1.5-2.5. It should be noted that different shading patterns of green, blue 
and gray are a consequence of overlapping shading and do not indicate additional information within the figure. The Galactic absorption column is indicated in 
the left panel by the horizontal dashed line. 
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FIG. 12. — Lx plotted against the flux ratio, where the 0.3—8.0 keV flux from the input two-component PO+DBB model (black line) is compared to the Lx 
values derived from the single-component fits; the PO model is indicated by the dashed-dot green line and the DBB model by the dashed blue line. Flux derived 
from a canonical PO model of Galactic JVh and T=1.7 is indicated by the dotted orange line. The panels on the left indicate the simulations with an input T=1.7 
and the right, the T = 2.1 simulations. Input AT; n values are labeled in each panel, as is the input photon value. 
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FIG. 14. — Best-fit results of simulations using three different photoelectric absorption models (detailed in section FOl following the same presentation as in the 
left-hand panel in Figure|9] Here, only the best-fit values derived from the T=1.7 simulations are presented for values of CTj n =0.50, 1.00 and 2.25 keV, over flux 
ratios values of 40% to 100%. The results from the phabs absorption simulations are labeled by the black diagonal lines, the results from the tbabs simulations 
are shown by the green lines and the results from the wabs simulations are plotted in blue. Error bars in the top left indicate the standard deviation for each set of 
simulations, these are also color-coded. The horizontal dashed line indicates Galactic Nn. 
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FIG. 15. — Best-fit results of simulations using four different Galactic absorption values (see section l4~4l fore more details). In this figure, instead of presenting 
the best-fit value of Nn, as shown in Figure|9] the value of excess A?h is plotted against the disc-to-total flux ratio, where the excess Nn is defined as the best-fit 
Nn value normalized by the Galactic absorption used in the simulations (Nu- Who). Here, only best-fit values derived from the T=1.7 simulations are presented 
for values of CTi„=0.50, 1.00 and 2.25 keV, over flux ratios values of 40% to 100%. The results from the main simulations are labeled by the black diagonal 
lines (with ^0=2.78 X 10 20 cm~ 2 ), the results from the NGC 4278 simulations (7Vho=1-76x 10 20 citT 2 ) are shown by the orange lines. The results from the 
/Vho=5 X 10 20 cm~ 2 simulations are plotted in blue and the JVho=1 x 10 21 cm~ 2 simulations are indicated by the green diagonal lines. Error bars in the top left 
indicate the standard deviation for each set of simulations, these are also color-coded. 
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FIG. 16. — A decision tree summarizing how the best-fit absorption values from the PO and DBB single-component models from the simulations can be used 
to determine spectral source states. In this figure the interpretation of each combination of /Vh(po) an d n H(dbb) is presented in the ellipses, where TD, NT, and 
intermediate source states are given, including sources with cool-disc components. 
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FIG . 17. — A population plot presenting the Lx and kT m or V values of the eight sources presented in this paper, where multiple points are plotted for spectrally 
variable sources to indicate the best-fit values from each spectral grouping. In this figure the left panel presents Lx~kT m for sources in a thermally dominant state 
and the right panel presents Lx— r for sources in a hard state. In both panels the points represent the results from the spectral fits, with 1<j errors indicated. Each 
source is denoted by symbol and color, with labeling provided in the bottom right corner of each panel. In cases where a source has been determined to be in 
an intermediate state dotted error bars are presented. In the case of observations 2&3 for S102 the thermal and non-thermal components are plotted separately 
in both panels and are indicated by dashed error bars. In the L\-kT m panel this point is not included in the main figure but within an insert, where, due to the 
cool temperature of the disc (0.14 keV), the x-axis has been adjusted to cover Ki m values between 0.1—1.0 keV. In the main L— kT m plot the diagonal dashed 
lines indicate the L oc7" 4 relation for BH masses of 5, 10, and 20 Mq , following the relation of Gierliriski & Done (2004). In the insert plot, in addition to these 
relations a line indicating L ocT 4 for a BH mass of 500 Mq is also presented. 



